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We propose an optical method for increasing the number 
of atoms in a pair of dilute gas Bose-Einstein condensates. 
The method uses laser-driven Raman transitions which scat- 
ter atoms between the condensate and non-condensate atom 
fractions. For a range of condensate phase differences there is 
destructive quantum interference of the amplitudes for scat- 
tering atoms out of the condensates. Because the total atom 
scattering rate into the condensates is unaffected the conden- 
sates grow. This mechanism is analogous to that responsible 
for optical lasing without inversion. Growth using macro- 
scopic quantum interference may find application as a pump 
for an atom laser. 
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In the recent experiments demonstrating Bose-Einstein 
condensation of alkali vapors the first stages of cooling 
are optical ||] . The final stage utilises evaporation of the 
hottest atoms out of the trap ||^ . Despite the great suc- 
cess of evaporation it has the disadvantage of removing 
atoms from the system. Consequently alternative final 
stage cooling methods are being investigated. Velocity 
selective coherent population trapping (VSCPT) is one 
optical method potentially capable of cooling to the Bose- 
Einstein transition point ^ . 

In this paper we analyze an optical method for increas- 
ing the number of atoms in a pair of overlapping Bose- 
Einstein Condensates (BECs), such as have been pro- 
duced in the laboratory using sympathetic cooling [§,1). 
According to conventional spontaneous symmetry break- 
ing arguments BECs are in coherent states with a def- 
inite global phase A pair of BECs therefore has a 
definite phase difference, which can be measured by a va- 
riety of techniques [[7|-pT[ . We utilise this phase difference 
to suppress transitions of atoms out of the condensates 
while driving transitions into the condensates, produc- 
ing a net condensate growth. The transitions are driven 
by spontaneous Raman scattering of laser light ^ . For 
a certain range of condensate phase differences transi- 
tions out of the condensates are suppressed by destruc- 
tive quantum interference. The interference is between 
the amplitudes for transitions from each condensate into 
the non-condensate atom fraction. 

This suppression of atom scattering out of the con- 
densates is analogous to the suppression of absorption 
that occurs in optical lasing without inversion p^ . In 
both cases the suppression is due to destructive quan- 



tum interference between two channels to the same final 
state. This mechanism for condensate growth might also 
be suitable as an optically driven pump for an atom laser 

Spontaneous Raman scattering of laser light from two 
BECs has been analyzed by Ruostekoski and Walls 
whose analysis we shall follow. They showed that the 
scattered light spectrum depends on the phase difference 
between the two condensates. The spectrum has two 
peaks, corresponding to transitions of atoms into and 
out of the condensates. For certain values of the con- 
densate phase difference the second peak disappears due 
to the destructive quantum interference previously de- 
scribed. Although this suggests the possibility of con- 
densate growth, only the short time behaviour was con- 
sidered. In fact, the Raman lasers induce Josephson os- 
cillations between the condensates [Q. We show that 
growth occurs and can persist over a complete period of 
this dynamics. However the necessary destructive inter- 
ference only occurs for a particular range of condensate 
phase differences. Consequently, condensate growth only 
occurs for a subensemble of condensate pairs. Suitable 
subensemble members might be chosen after a measure- 
ment of the phase difference 0-0. 

We consider the set-up first introduced by Javanainen 
0; a quantum degenerate gas with two ground states 
and a common excited state. All the ground state atoms 
are confined in the same trap. The pair of condensates 
are in two different Zeeman sublevels \b) = \g,'rn) and 
|c) = \g,m — 2). The state |c) is optically coupled to the 
electronically excited state |e) = |e,m — 1) by the field 
E2 having a polarization a+ and frequency Similarly, 
the state \h) is coupled to |e) by the field Ei with a po- 
larization (T„ and frequency Vli. Following Ref. [|l^ the 
Hamiltonian density for the system is 

- (db ■ El VlV'e + d, • E2 ibl^Je + h.C.) . (1) 

The first three terms reflect the center of mass energy, 
H, and the internal energies of the atoms in the absence 
of electromagnetic fields. The frequencies for the optical 
transitions e •«-> 6 and e c are tUet and Wgc i^cb = ^eb — 
ujec)i respectively. Ti^ is the Hamiltonian density for the 
free electromagnetic field. The final, bracketed, terms are 
for the atom-light dipole interaction. The dipole matrix 
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element for the atomic transition e <-> 6 (e <-> c) is given 
by df, (dc). 

We assume that the driving Ught fields are in coher- 
ent states and detuned far from single photon resonance 
so that multiple scattering can be ignored. We also as- 
sume they are plane waves propagating in the positive z 
direction with wavevectors Ki 



Ed2(i") = -^^26+ exp(iK2 • r) , 



(2a) 
(2b) 



bch = ^\—^2~ ^cb^ and the Raman coupling coefficient 

K 



Si = 
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The dipole matrix elements df, and dc contain the reduced 
dipole matrix elements and the corresponding nonvanish- 
ing Clebsch-Gordan coefficients. To simplify the algebra, 
we assume k to be real. 

The intensity of the scattered light at position r is 
given by 



where are the unit circular polarization vectors, and 
we have defined slowly varying fields by = e'^'*E^. 
We also define the slowly varying matter field ipc — 

In the limit of large detuning the excited state field 
operator may be eliminated adiabatically. Following 
Javanainen and Ruostekoski the scattered electric 
fields may then be expressed in terms of the driving fields 



/(r) = 2ceo(E;.E+), 
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E^2 where E is radiated by decays into 



X jdb ■ E+ ^lijb + d, • E+ Vl^c} 



(3) 



The driving fields and atom fields are all functions of 
r' and t. Ai = ili — ujeb is the atom-field detuning of 
field 1. The first line represents the radiation from the 
atomic dipole density, and the second follows after adi- 
abatic elimination of the excited state. E^2j which is 
radiated by decays into state |c), is found by swapping 
subscripts b and c and swapping the driving fields EJ-^ 
and E^2- We have used the first Born approximation 
based on the assumption that the incoming fields domi- 
nate inside the sample, as multiple scattering is negligi- 
ble. The kernel K(d) is the familiar expression |16| for 
the positive-frequency component of the electric field at 
r from a monochromatic dipole with the complex ampli- 
tude d, located at r'. 

After adiabatic elimination of the excited state from 
the Hamiltonian density Eq. (|l|), and approximation of 
the electric fields by the driving fields, the following 
Hamiltonian density is found to first order in the inverse 
atom-field detuning |^ 

Hm = ^1{H ~ Mi)^f, + ^1{H ~ hS,i, - h52)i>c 

+hK [i^l-ipc exp(-iKi2 • r) -I- ■0j'06 exp(i/«i2 • r)^ , (4) 

where K12 = Ki — K2 is the wavevector difference of the 
driving light fields. We have introduced the light-induced 
level shifts 5i, the detuning from two-photon resonance 



Substituting in the expressions for the scattered fields 
in terms of the atom fields, Eq. (^), generates a sum of 
terms for the intensity of the form 



2ceo (^^) j dV'dV"[K(db)r ■ K(df,)" 

x(d^ E,- )(dfc . E+ )(vi Vv-rV) • 
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The ' and " respectively denote functional dependence 
on r' and r". We now assume that the driving fields 
have the same wavevectors, so that H12 = 0. The 
dynamics of the ground state fields tpc and tpb follows 
from the Hamiltonian Eq. (^) . We assume a translation- 
ally invariant and non-interacting Bose gas. The mat- 
ter field operators are given by the familiar plane wave 
representations i/jbi'rt) — V~^^^ ^^exp{i'k ■ r)6k(i) and 
'ip(.{rt) = V~^l'^ exp(ik • r) Ck(t), where V is the mode 
volume. In the absence of light, the center of mass mo- 
tion in both ground states satisfies the dispersion relation 
Ek = ?i|kp/2m, with m the atomic mass. Defining the 
effective two-photon detuning 28 = bch — §\ + ^2 and the 
condensate oscillation frequency fi^ = (S^ + k^)^/^, the 
mode operators at time t are given in terms of the oper- 
ators at time t = by 



5k (t) = e"* {A£k(0) - B6k(0)} , 
6k(t) = e'"* {A*fok(0) - i?5k(0)} , 



a = (5 + (5i 
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A — cos VlYit + — — sin B 



(8a) 
(8b) 
(8c) 



i-——smiljit. (8d) 



Before the light is switched on, the atoms in the states 
\b) and |c) are assumed to be uncorrelated. The Raman 
fields induce a coupling between the two levels. In the 
presence of Bose condensates in the ground states, the 
coupling between the two condensates is analogous to 
the coherent tunneling of Cooper pairs in a Josephson 
junction 0,0. The expectation values of products of 
four atom field operators, such as occurs in Eq. (0), mav 
be evaluated after substituting in the expressions Eq. (g) 
for the fields at time t in terms of the time zero fields. 
For example the expectation value in Eq. (|^) becomes 
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where all the field operators on the right hand side are 
evaluated at time zero. 

The field operators are sums over the condensate and 
non-condensate modes. Since we are only interested in 
the change in the number of condensate atoms due to 
light scattering we need only evaluate those terms corre- 
sponding to scattering of atoms into or out of the conden- 
sate, i.e. the incoherent part of the scattering. We ignore 
scattering of atoms between non-condensate modes. To- 
gether with momentum conservation this leads to a con- 
siderable simplification of the terms like Eq. (^. Once 
a particular plane wave mode is chosen for the first fac- 
tor in Eq. (|9|) the requirement for a non-zero expectation 
value determines the modes occurring in all the remain- 
ing factors. For example, the part of Eq. (^) relevant to 
condensate depletion and growth is 



{dId^dIdo 
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where Di — A*{t)b,i,{0) — B{t)ci,{0) and the subscripts 
and ± respectively refer to the condensate mode and 
the non-condensate modes having momenta ±hAK. Here 
Ak — rin/c—K is the wavevector change of the scattered 
photon, and n = r/|r| is the unit vector in the light 
scattering direction under consideration. The two non- 
condensate modes +/— respectively arise from scattering 
of atoms into/out of the condensate. Note that the par- 
ticular atomic mode denoted depends on the light scat- 
tering direction n, as does the polarization of the scat- 
tered light. In general, the polarizations of the emitted 
photons from the two different atomic transitions are not 
orthogonal. However, although the resulting interference 
terms are nonvanishing in a particular direction ^ , their 
contribution to the total intensity of the scattered light 
vanishes after integration over all scattering directions. 
Because we are ultimately interested in the total inten- 
sity, we ignore the terms proportional to (E^^ -E^j) +c.c. 
in the intensity of the scattered light. 

For brevity we assume that the number of atoms in 
the ground non-condensate states arc the same, n± — 
{bj.b±) — (c]^c±), and that rt+ = n_ = n due to isotropy. 
Further simplification occurs if we assume that there are 
equal numbers of atoms N in each condensate, and that 
the laser intensities are chosen so that the level shifts are 
equal S = Si — $2. Evaluating all the relevant terms in 
Eq. ^ we find the following expressions for the intensity 
due to scattering of atoms into and out of the condensates 



/in = 2CnS(f,{N + 2Re[A*B{clbo)]dl} , 
/out = 2Cin + l)Sdl{N + Re[{A*^ - B^){clbo)]} . 



C = 



STT^eoAi |r| 
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(11) 
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(13) 
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and 



where Cl is the speed of light, d^ 

{d^ — d^)/(/^. We next integrate these intensities over the 
Josephson oscillation period P — 2-K/VtB.. We find the 
time averaged intensities 
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C'n <! 1 + -^d^cose 



/outrfi = C"(n + l)|l + ^cose| 



(14) 
(15) 



where C = 2C6dlN, and 9 is the condensate phase dif- 
ference. The non-condensate populations n are functions 
of Ak and hence of the scattering angle 9 between n and 
the laser propagation direction. Our final step is integra- 
tion over all scattering directions. This yields the total 
scattered light intensity and hence the total atom tran- 
sition rates. The angular integration has the effect of 
replacing n by 



Q /-TT -1 

— n = 27r / {l--sm'^9)n{9)sinede, 
3 Jo 2 



(16) 



and n -|- 1 by 87r(Ti + l)/3. This integral may be inter- 
preted as the number of non-condensate atoms available 
for scattering into the condensates. In an infinite homo- 
geneous system the integral is divergent at the low energy 
end. However, for a finite, trapped system a low energy 
cutoff is provided by the first excited state. A numer- 
ical integration assuming the Bose-Einstein distribution 
at T = 400 nK, and a low energy cutoff of ?i(27r x 100) 
J, gives (87r/3)n w 65 for rubidium. However, this is a 
crude estimate since realistic systems are not expected to 
be in thermal equilibrium. 

Assuming, for simplicity, equal dipole moments dd = 
the net rate of scattering of atoms into the condensates 
(atoms per second) is then, from Eqs. (|^-|l^). 



R = 



Id 



N 



K 



1 TT2- COS e [1 



(17) 



where 7 = dl\K\'^ / {Sneoh) is the free space spontaneous 
emission rate, and Id is the intensity of the lasers. Con- 
densate growth corresponds to a positive rate R. As- 
suming that n 3> 1 this is equivalent to the following 
requirement on the condensate phase difference 



1 

cosO < --^ . 
n K 



(18) 



This inequality is fulfilled by particular negative values of 
cos 6 provided that n is sufficiently large, and that the ef- 
fective two-photon detuning 26 is sufficiently small. The 
latter may be chosen small by manipulating the light- 
induced level shifts or the relative frequency of the driv- 
ing light beams. The assumption of equal wave numbers 
for the driving light fields in the calculations is not very 
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restrictive for the relative frequency, because an atom 
trap introduces an uncertainty for the momentum con- 
servation. In the limit of small two-photon detuning the 
effective linewidth of the transition c — > 6 may have an 
effect. However, it may be shown to be proportional to 
Aj^^ or smaller. 

With a 1 /im wavelength, 1 mW cm^^ laser intensity, 
and (K2/0|j)cose = -1/2, the growth rate Eq. (f?!) is 



R^Nn 



(19) 



A laser detuning of A1/7 = 10'^, condensates with 
N — 10^, and {8TT/3)h = 65 give a condensate growth 
rate of i? w 7 x 10^ atoms per second. This rate is large 
enough to be useful for both atom laser pumping and for 
condensate growth. However, sustained growth will re- 
quire repopulation of the relevant non-condensate atom 
modes by atom scattering processes. Other limitations 
on growth include diffusion of the condensate phase dif- 
ference due to atom-atom interactions Jl7| , p^ , and heat- 
ing of the non-condensate atom fraction by Raman tran- 
sitions. It is also required that the spatial overlap of the 
condensate wave functions is significant. 

Multiple scattering can be ignored provided that only 
a small fraction of the incident photons are scattered. If 
the cross sectional area of the condensate is A, then this 
is true, if > 67r(7/Ai)2iV7i/|/«p. With the preceding 
parameters this becomes ^ ^ 4 x 10^^/xm^, which is 
easily fulfilled. 

In the calculations we used the conventional sponta- 
neous symmetry breaking arguments that BECs are in 
coherent states. However, this convenient approach is 
by no means necessary. In fact, for the present atomic 
level scheme the relative phase between the two conden- 
sates has been established in stochastic simulations of 
the measurements of spontaneously scattered photons, 
even though the condensates are initially in pure number 
states . Without any measurements of the condensate 
phase difference the macroscopic quantum coherence is 
expected to undergo collapses and revivals (l8| . Because 
the detections of spontaneously scattered photons estab- 
lish the relative phase, they could also stabilize the phase 
against the collapse of the macroscopic wave function. 

We have shown that quantum mechanical interference 
enables growth of two Raman driven Bose condensates. 
The growth mechanism is analogous to that for gain 
in the laser without inversion. In that case destructive 
quantum mechanical interference suppresses absorption, 
allowing gain to dominate. In our case interference sup- 
presses scattering of atoms out of the condensates, allow- 
ing scattering of atoms into the condensates to dominate. 
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